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We study the impact of a generalized Chaplygin gas as a candidate for dark energy on CMB 
anisotropies. The generalized Chaplygin gas is a fluid component with an exotic equation of state, 
p = —Ajp'^ (a polytropic gas with negative constant and exponent). Such component interpolates 
in time between dust and a cosmological constant, with an intermediate behaviour as p = ap. 
Perturbations of this fluid are stable on small scales, but behave in a very different way with respect 
to standard quintessence. Moreover, a generalized Chaplygin gas could also represent an archetypal 
example of a phenomenological unified models of dark energy and dark matter. The results presented 
here show how the CMB anisotropies induced by this class of models differ from a ACDM model. 


PACS numbers: 98.80Cq 

1 . According to the present recipe to explain obser¬ 
vational data, two dark components seem to fill the uni¬ 
verse up to 95% of the total content. The baryon density 
is indeed very small, flbh'^ ~ 0.02 |^, assuming a flat 
Universe. This situation has been slowly reached over 
decades. In addition to cold dark matter (CDM), in the 
90s a cosmological constant term A was called into play 
to explain the recent accelleration of the universe indi¬ 
cated by the Supernova data |Q, but then confirmed by 
other observations. 

Avoiding to explain on theoretical grounds the em- 
barassing smallness of a cosmological constant A con¬ 
strained by observations, a scalar field Q [||, dubbed 
quintessence, was su^ested in order to explain the ac- 
cellerating Universe [Q. Because of the different back¬ 
ground evolution and of the presence of fluctuations, 
the challenge is to distinguish A from quintessence, for 
instance through the Cosmic Microwave Background 
(CMB) anisotropies |^. 

An alternative to quintessence for modelling a dark 
energy component could be a perfect fluid (henceforth 
PF) with a generic pressure p = p{p) (not being linear in 
energy density p as for barotropic fluids) whose energy 
momentum tensor is: 

= Pg^LV + [P+ P)Uf_,Uy ( 1 ) 

where g^^ is the metric and is the fluid velocity 
{Ufj,u^ = —1). From a theoretical point of view a scalar 
held description would be preferrable in the logic to re¬ 
late an accelerating universe with a fundamental quan¬ 
tum. From a phenomenological point of view the reasons 
to prefer one over the other are less obvious. An ex¬ 
otic fluid capable to develop a negative pressure at late 
times may represent as well the effective degree of free¬ 
dom which drives the present accelleration of the uni¬ 
verse. In particular, a PF model with this property would 
allow to explore the possibility that dark energy clusters 
on small scales. Indeed, it is important to note that, 
by parametrizing dark energy with an uncoupled scalar 


field with ordinary kinetic term, one does implicitly as¬ 
sume no clustering of dark energy on scales smaller than 
the Hubble radius. 

Among the class of PF models which could work as a 
dark energy component - in principle one could design 
suitable pressure profiles p{p) instead of potentials V ((/>) 
for a quintessence field (f - the Chaplygin gas has 
recently received a lot of attention 0. A Chaplygin gas 
is characterized by a pressure px related to the energy 
density px in the following way: 


Px = — 


A 

Px 


( 2 ) 


with a = 1. Even if such exotic fluid was proposed in 
the context of aereodynamics there are interesting 
connections with particle physics and d-branes ||]. A 
Chaplygin gas is also equivalent to a tachyon field with a 
constant potential [|| and, at the homogeneous level, to 
a complex scalar field |j^ or to a quantum scalar held 
|]lT| in the bottom of a potential (called Thomas-Fermi 
approximation in j^). In this paper, we study a gen¬ 
eralized version of the Chaplygin gas (henceforth GCG) 
I p^ by considering 0 < a < 1 in Eq. (||). 

2 . In a Robertson-Walker metric 

ds^ = -dt^ + -V r'^dQ'^) , (3) 

1 — Kr^ 


where K = 0, ±1 is the curvature of the spatial sections 
and 0 is the solid angle, the energy conservation equation 
for a GCG 


px + iH {px +px) = 0 


can be immediately integrated [[7|l2|: 


PX 



B 

a3(i+«) 



( 4 ) 

( 5 ) 


where A, B are constants with dimensions We 

note that a GCG reduces to a ACDM model for a = 0 
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and to a sCDM for A = 0. We see that this fluid with 
an exotic equation of state behaves like dust for small 
a (when B/A » assuming a = 1 at present 

time) and a cosmological constant given by in the 

opposite limit {B/A << By Taylor expanding 

in this limit we obtain from Eqs. (||) and (^): 


(^+ (l + a)Aa3(i+a) +^(^ 12 )) 

^ ^ (1^)) 


( 6 ) 



-10 -8 -6 -4 -2 0 


Therefore, in this limit a GCG behaves as a sum of a 
cosmological constant and a perfect fluid characterized 
by p = ap, sheding light upon the physical meaning of 
the parameter a. We note that the solution (||) to Eq. (^) 
is valid for any a > —1, but also for a < —1 (a standard 
polytropic gas). In this latter case the behaviour of such 
a PE interpolates between a cosmological constant and 
dust. Eor a = — 1 Eq. (||) describes the usual barotropic 
perfect fluid. According to Eq. (^) , the equation of state 
wx defined as: 


Px A 

^^ = — = — 
Px Px 


(7) 


decreases from the value 0 to —1. An example of a back¬ 
ground evolution for a GCG as dark energy is given in 

Fig. ( 0 ). 

The presence of a GCG as a dark energy component 
could be distinguishable from a quintessence component 
because of the parametric form of the pressure p = p{p). 
The time derivative of the equation of state of the dark 
energy component is important in the program of recon¬ 
structing the total equation of state m- Indeed, the 
time derivative of wx for a GCG is 


Wx = 5H{a + l)wx{l + Wx) ( 8 ) 

while for a scalar held (/> with potential V = V (</>) is: 


= 3H{1 + w^){w^ - 1) - 2—. (9) 

P4> 

SN la observations can constraint the GCG as a can¬ 
didate for dark energy, as recently studied by different 
authors |^]. The purpose of this letter is to show how 
CMB data could be more selective. 

3 . The behaviour of perturbations is a very interest¬ 
ing aspect of the model. When trying to build a model 
for an accelerating universe with a barotropic (constant 
equation of state w) perfect huid, one runs in the prob¬ 
lem of instabilities on short scales because of a negative 
sound speed for perturbations. Infact, the sound speed is 
equal to the equation of state w and this must be nega¬ 
tive {w < —1/3) to explain the accelleration. This is the 
usual problem for a fluid description of domain walls and 
cosmic strings. Quintessence models with scalar helds 


FIG. 1. Evolution of the background quantities Q.r (yel¬ 
low), (red), Qx (blue), Wtot (purple) and wx (green) for a 
flat universe versus the scale factor. The horizontal black line 
is —1/3, which denotes the threshold for w below which the 
universe expands accellerating. The parameters are h = .7, 
Qro = 10■^ Qmo ^ .3, Qxo = .7, Q = 1, B/A = .01. 


with a standard kinetic term do not have this problem. 
The sound speed for a scalar held is equal to 1. K-essence 
models based on scalar helds with a non-standard ki¬ 
netic terms are different in this respect , but still have 
positive sound speed. For a GGG the sound speed Cx for 
perturbations is 

2 9p A 

^x = -^ = an+T = ■ (10) 

dp 

Therefore, because of its non-barotropic nature, pertur¬ 
bations of a GCG are stable on small scales even in an 
accellerating phase, and behave similarly to dust pertur¬ 
bations when the gas is in the dust regime. When the 
behaviour of the background Chaplygin gas is of A-type, 
the sound-speed is a. In order to avoid causality issues, 
we shall consider a < 1. This latter constraint, with 
the requirement of positive sound speed (a > 0), marks 
the interesting physical range of a. We note that the 
possibility of having a non negative sound speed and a 
negative equation of state, as happens for a GCG, could 
open a new development in modelling a domain wall or 
cosmic string network which is not plagued by short scale 
instability for perturbations. 

The equations for the energy density contrast Sx = 
5px /Px and the velocity potential Ox in the synchronous 
gauge are, according to Ref. and by using Eqs. ([sHl^) 

5'x = —(1 + Wx) (^Ox + —^ + 3H{wx — Cx)Sx 

0'x = -n{l-3cj,)ex + -^^k^Sx, ( 11 ) 

1 -I- wx 

where h is the trace of the metric perturbations in the 
synchronous gauge |j^ . This set of equations agrees with 
those used in ||^. In order to study the Jeans instability 
for a GCG it is useful to study the equation for the (gauge 
invariant) comoving density contrast Ax = Jx +3(1-1- 
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FIG. 2. Density contrasts of photons (solid line), baryons 
(dashed line), GCG (short dashed, on the left) and GDM 
(short dashed, on the right) in a GCG model (left) and a 
ACDM model (right). 


wx)'H9x/k‘^ in the approximation in which the GCG is 
the only component of the universe: 

Ajjf +7^(1 + 3c^ — 6tcx) 

+ cxk'^Ax - (l + Swx - - 6cx) Ax = 0. 

We see that the Jeans instability of a GCG is very similar 
to GDM in the dust limit (when wx ^ c\ ~ 0). How¬ 
ever, because of the time dependence of wx the Jeans 
instability is progressively removed since the quantity 
1 -I- Swx — Swx — Gcx changes sign as wx departs from 0. 
At the same time, GCG perturbations start to oscillate as 
soon as draws away from 0. We confirm numerically 
this behaviour, also when other fluids are present. In Fig. 
(^ we show the comparison of the evolution of cosmolog¬ 
ical perturbations in GCG models (without CDM) and 
ACDM models. 

We have implemented Eqs. and ( 0 ) in a modi¬ 
fied version of CMBFAST Q . We have tested the code 
against the sCDM model obtained by setting A = 0 in 
Eq. ® and we have considered an initial adiabatic scale 
invariant spectrum for perturbations. Instead, by set¬ 
ting i? = 0 in Eq. (||) one obtains a ACDM model in 
the background, but with not-trivial perturbations in the 
dark energy sector (see Eq. ®)- 

In Fig. (^ we show the dependence of the spectrum of 
temperature anisotropies on the ratio B/A and a. The 
three parameters A, H, a are in direct relation with the 
physical quantities of dark energy at present time, as 
wxQ,wxo- Therefore these parameters can be con¬ 
strained by maximizing the likelihood function with the 
present data |^. In particular, from the left panel of 
Fig. i) one can see how spectra can be sensibly differ¬ 
ent for a GCG model which differs from a ACDM model 
by less than 10% in the background evolution. 

Because of its early dust behaviour, a GCG may also 
represent a prototypical unified model of dark matter and 
dark energy (see |^| for a similar proposal, but with a 



FIG. 3. Ci spectra of temperature anisotropies 5TjT ver¬ 
sus I for GGG varying the ratio B/A (on the left) and a (on 
the right). The parameters are the same of Fig. 1 except 
B/A = 0.1,0.01,0.001 and AGDM (from bottom to up, on 
the left) and a = 1, 0.5, 0 (from bottom to up, on the right) 
for B/A = 0.1. 



FIG. 4. Spectra of temperature anisotropies 5T for a uni¬ 
fied model (Djifo = -95 and h = 70) varying B/A compared 
with the BOOMERANG, MAXIMA, DASI data. The solid 
line at the bottom is B/A = 1. The blue lines correspond 
to 10, 50 and a sGDM respectively from bottom to top. The 
red lines correspond to .01, .001 and a A-baryonmodel from 
bottom to up. 


scale dependent equation of state). In Figs. (HJI) we 
compare GCG models without CDM {h = .7 i ffxo = -95) 
with the BOOMERANG ||, MAXIMA || and DASI 

f data, when varying B/A and a. In particular, in Fig 
we see how all the models lie below the limiting case 
of sCDM (A = 0) and a A-baryon model (very similar, 
but not equivalent to H = 0). The resulting spectrum 
of CMB anisotropies constrain these GCG models |2^] 
much more than the SN la data . 

4. We have studied the implications on the evolutions 
of cosmological perturbations and on CMB anisotropies 
of a GCG as a candidate for dark energy. This GCG 
covers all the interesting possible cases of a dark energy 
model from a polytropic gas. A GCG is more distin¬ 
guishable from a ACDM model than a QCDM model ||] 
since both GCG background and perturbations are im- 
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